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,__! i Abstract. We study the convergence of the series expansions sometimes used in the 

CO I analysis of the nuclear effects in Deep Inelastic Scattering (DIS) proccesses induced 

' by leptons. The recent advances in statistics and quality of the data, in particular 

f~i . for neutrinos calls for a good control of the theoretical uncertainties of the models 

I I used in the analysis. Using realistic nuclear spectral functions which include nucleon 

^ • correlations, we find that the convergence of the derivative expansions to the full results 

3 , is poor except at very low values of x. 
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O ■ 1- Introduction 

vn ; 

^ I The advent of new and high statistics deeply inelastic neutrino scattering experiments 

00 ! apart from providing valuable data on F2 and a; F3 [H [2] has shown again the importance 

T\ i of nuclear effects that render difficult the extraction of the nucleon structure functions. 

.^ \ Furthermore, the use of nuclear targets is necessary due to the smallness of the cross 

rS ' sections. Qualitatively, the nuclear effects are well known. Shadowing, antishadowing, 

c^ I Fermi motion, binding,... have all been widely studied for charged leptons in the context 

of the European Muon Collaboration (EMC) effect. For a review, see [31 Hj. Moreover, 

it is well known that nuclear effects are substantially different for neutrino reactions due 

to the presence of the axial current and the different valence and sea quark contributions 

for each observable [5j. This situation asks for a detailed and quantitative microscopical 

understanding of the nuclear effects, rather than the parametrizations that have been 

used sometimes, like recently by the NuTeV Collaboration [2]. 

One of the basic ingredients in all calculations is the nuclear spectral function. 
This presents some serious difficulties, as these functions are not so well known to 
the precision level reached by current experiments. Therefore, the analysis might 
introduce unwanted model dependences. However, it was soon noticed that under some 
approximations the nuclear structure functions could be written as simple expansions 
on the nucleon structure functions and their derivatives. All the nuclear information 
would then be encoded in the expected values of some nuclear magnitudes, like the 
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average kinetic energy of the nucleons or the mean nucleon removal energy. These kind 
of approximations, both for charged lepton and for neutrino induced reactions, have 
been widely used in the literature P, [3, El El EDI [HI [121 US] • The aim of this paper is 
to investigate the quality of these expansions and to what extent they could be used in 
the analysis of lepton deep inelastic scattering experiments. In order to do that, once 
the formalism is established in the next section, we will make a comparative study for 
^2^3 and for a few typical nuclear spectral functions used in the literature. 

2. Formalism 

The nuclear structure functions can be written as a convolution of the nuclear spectral 
functions and nucleon structure functions. See for instance [Hj and references therein. 
In the rest frame of the nucleus, the F2 and F3 structure functions are 

T=p,n I \ ^ ^ / 

FiM') = E /^^^('.P) fl + ^) ^FU.',Q% (2) 

T=p,n ' 



(27r)4 ^ '^-^ \^ ^Mj X 



where Vp^"-' (e, p) is the nuclear spectral function, normalized to the number of protons 
(neutrons) in the nucleus, and describes the probability of finding a proton (neutron) 
with momentum p and removal energy e. The four-momentum of the nucleon can be 
written as p = (M + e, p), with e < 0. The z axis is oriented in such a way that q lies 
on it, p_L is the transverse momentum of the nucleon and 7 = |q|/5''^. Here, x' is the 
natural Bjorken variable for the nucleon in the nuclear medium, i.e. x' = Q'^/{2p ■ q); 
while X is the Bjorken variable in the nucleon rest frame, x = Q'^/{2Mq^). They are 
related by 

x' = — where z = 1 + -— + 7-7^- (3) 

z MM ^ ^ 

For isoscalar nuclei such as '^^Ca, only the isoscalar component of the spectral 

function and the structure function have to be accounted for. When this is done ([H) 

and ([2]), these read as 

..,..«^)^./gn(.p)^(l.4f!.6^).."(f.«^). (4) 



ded^p f p 



fJy , AT / tXj 



F.\^,Q^) -A J ^n(.,P) ^1 + :^j ^Ff [^,Q') , (5) 

where Po(e, p), which is the isoscalar part of the nuclear spectral function, is now 
normalized to unity and we perform the calculations for the nuclear structure functions 
averaged over neutrinos and antineutrinos, i.e., we only consider the symmetric u + u 
combination P^ [T5] . 

In a nucleus, the expected values of ^ and -j^ averaged with the nuclear spectral 
function are much smaller than unity. Thus, z ^ 1 and x' ~ x. Under these assumptions. 
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we can perform a Taylor expansion of the integrands in expressions (jl]) and ([S]) around 
z = 1, keeping terms up to order e/M and p^/M^. In this way we will be able to take 
out of the integral the structure functions and their derivatives and we will be left with 
expected values of the removal energy e and momentum squared p^. This statement 
is true if the nuclear spectral functions only depend on the modulus of the momentum 
IpI (as it is in the case of the nuclear spectral functions we will consider) and not upon 
its direction. Under this assumption we can drop the expected values of p^ (or any 
other component of the momentum) and epz because they are identically zero due to 
symmetry considerations. 

We will begin with the structure function F^ (x, Q"^). After performing the Taylor 
expansion, keeping terms up to order e/M and p^/M^, and dropping those terms which 
go with {pz) /M or {epz) /M^ (because they are identically zero as stated above), we 
obtain 



Fi 



a;,Q2 



A 



F^{x,Q') 



l+(7^-l) 



2\ H 



+ 



X 



dF^ 



x.Q' 



+ 



x^ d^Ff 



M 

3M2^ 
x,Q2 



— X- 



dFi' (x, Q' 



1 



dx 

(p') 



M ^ ^ ' 3M2 



(6) 



dx 2 dx^ J ' 3M2 

where (O) is the expected value of the operator 0{e, p) averaged with the isoscalar 
spectral function Vo{e, p): 

ded^p 



{O) 



{27ry 



Po(e,p)0(e,p). 



(7) 



We have also used the fact that, with a nuclear spectral function which depends only 
on the modulus of the momentum, the expected value of momentum squared is shared 



equally among every squared component, i.e: (pf) 



(P^) 



In the Bjorken limit, (7 



m- 



-^3 



coincides with the expansions used by Frankfurt et al. [6] and Ciofi degli Atti et al. 

If we do the same for the nuclear structure function F^ (x, Q"^), we obtain 
(^>Q')-,^7v/ ^2^ {^)\r.Nt ^2^ , dFf{x,Q') 



A 



+ 






X 



3M2 



1) 



Ff + x 



dFf 
dx 



+ 7 



X 



dx 
dx 



x^ d^F. 



N 



+ T 



dx"^ 



(8) 



a similar result to that Kulagin [10]. To allow for an easier comparison, in the above 
expressions, 7^ can be rewritten as 



7 



1 + 



AM^x 



2^2 



Q' 



(9) 



Then, the only difference between ([S]) and that of [TU] is that, for simplicity, we do 
not consider the off-shell dependence in the nucleon structure function F^ . In our 
calculation, the only source of off-shell dependence is through the nuclear spectral 
function Vq. Thus, we do not have the term dp2Ff that appears in [TU] . 

For our study we have selected the ^^Ca nucleus, which is isoscalar and already will 
show important medium effects. The free nucleon structure functions have been taken 



Study of the derivative expansions for the nuclear structure functions 



X 



^ 



i.O 
















X 




















1 ' 










(I) Full 














- 


- 


- 


(II) Full 


























I 


1.2 




■ 


■ 




(III) Full 


























n 


^ 


































/ T" 






































If 






































If 






































If 
if 


1.1 


- 


































/ jf 


1 


=1-3; 


^ 




































_ 










"""■■ 


"— 


_L 


■=■ 


— 








— 


— 


r-r 


■r-i 




1 


0.9 













0.2 



0.4 



0.6 



X 



Figure 1. R2 = ^ ratio for ""^Ca at Q^ = 20 GcV^ with tlic nuclear spectral 
functions (I), (II) and (III) described in the text. 



from [T6] and we have chosen several different nuclear spectral functions. The first one 
(labelled I) is a phenomenological model of the spectral function which has a mean field 
part and high-momentum components coming from NN-correlations and it is described 
in Kulagin et al. [^ where it was used in a global study of nuclear structure functions. 
The second one (labelled II), that also contains correlations is taken from [17] where it 
has been tested in the calculation of several electron scattering observables. The third 
spectral function is taken from the semiphenomenological model based on a many body 
calculation and the local density approximation that is described in [18j (labelled III). 
Finally, we also consider the simple mean field spectral function of [I2] (labelled IV), 
which was used in the study of the EMC effect. 



3. Results 



In figure [Tj we show the results of the ratio R2 = -^ for the nuclear spectral functions 
(I), (II) and (III). This gives us some idea of the uncertainties related to these functions. 
The differences are small, even when the spectral funtions have been obtained with 
diverse methods and are in fact quite different if one studies in detail their energy and 
momentum dependence. However, the expected values of the mean removal energy (e) 
and the mean kinetic energy per nucleon (T) = ( ^ j are quite similar, as can be seen 
in Table [H In particular, spectral functions (I) and (III) that have quite close expected 
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Table 1. Expected values of the nuclcon removal and kinetic energies for the nuclear 
spectral functions of Kulagin et al. [1^1 (I), Ankowski et al. [11] (II) and Fernandez de 
Cordoba et al. [THl (III). 



Spectral Function 


(I) 


(11) 


(III) 


(e) (MeV) 


-49 


-40 


-47 


(T) (MeV) 


30.1 


26.2 
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Figure 2. i?2 = ^ ratio for ^°Ca at Q^ ^ 20 GeV^. Comparison of the full results 
with the approximation of ([6]) for the nuclear spectral functions (I), (II) and (III). 



values also produce very similar ratios, whereas (11) which has an appreciably smaller 
binding energy gives slightly larger values for the ratio. 

In figure Ej we compare the full results for the same ratio (jlj) whith those obtained 
making use of the approximation of ([6]). The use of this expansion has been assumed 
to be a good approximation for x < 0.5 [ZIEIE]- Obviously, the series expansion agrees 
well at low values of x with the full results. However, in the three cases produces lower 
values for the ratio at intermediate x showing a maximum deviation of a 3-4% around 
X = 0.5 — 0.6. This region, with a dip in the ratio, is dominated by the mean removal 
energy (or equivalently the binding energy) per nucleon. Although this could look a 
small error, we should remark that it means increasing the deviation due to nuclear 
effects from the value 1 by around a 30% . At higher x's, where Fermi motion of the 
nucleons provides the dominant effect, the series expansions grow faster than the full 
results and become larger for x > 0.65. Thus, we find that for typical nuclear spectral 
functions the convergence of the series expansion is not so good except at very low x, 
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Figure 3. i?3 = ^ ratio for ^°Ca at Q^ = 20 GeV^ Comparison of the full results 
with the approximation of ^ for the nuclear spectral functions (I), (II) and (III). 



where in any case other effects not considered here, lilce shadowing play a major role. 

The results for the R^, = -^ ratio are shown in figure [3l The full model for R^ 
is given by ([5]) and the series expansion by ([8]). As it was the case for i?2, there is 
a dip region dominated by the binding energy and the Fermi motion of the nucleons 
produces the large rising at high values of x. The comparison of the full results and 
the series expansions shows the same features as for i?2- The expansions systematically 
underestimate the ratios at intermediate values of x and overestimate them for x > 0.7. 
This overestimation of the effect of the Fermi motion was already discussed in |10j . 
There, it was also claimed that in the limit of high Q^ and for heavy nuclei the expansion 
should be a good approxiamtion up to x < 0.75. However, after studying R2 and -R3 
at different Q^ values, apart from the one shown in this paper, we have found that for 
medium nuclei the convergence of the series is only good at very low x where other 
nuclear effects are very relevant. 

In order to obtain a better convergence one would need to reach a higher order in the 
expansions. This implies also the sensitivity to nuclear expected values of higher powers 
of the nucleons momenta as < (p/M)^ >. One should notice that due to correlations 
these expected values are not negligible [201 [19]. The three spectral functions considered 
above contain high momentum components and one may expect any expansion up to 
order p^ to fail to have a good convergence. To test this point, we have also calculated R2 
for the case of a simpler mean field spectral function which does not incorporate nucleon 
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Figure 4. i?2 = ^ ratio for "^^Ca 
approximation of 



-p2- icii-ii^ i"i ^^- Comparison of the full results with the 
for the nuclear spectral functions (IV). 



correlations and does not have those high momentum components [12]. This function 
was used to analyse the A dependence of the position and magnitude of the dip. Here, 
we have used the same parton distribution functions as in the original reference |12j . 
The results are shown in figure HI In this case, there is an almost perfect agreement 
between the full results and the expansion. 

In summary, we have studied the quality of some series expansions commonly used 
to incorporate approximately the nuclear effects in the analysis of DIS processes. We 
have found that for realistic enough nuclear spectral functions, that include nucleon 
correlations and have high momentum components, the convergence of the series to the 
full result is poor except at very low values of x, where in fact, other nuclear effects, 
like shadowing or antishadowing, are more relevant. At high x values, in the Fermi 
motion region, the expansions clearly overestimate the full result. This was known and 
expected. However, even at relatively low x values, where the expansions were suppossed 
to provide a good approximation, we have found that they systematically underestimate 
the value of F2 and F3, artificially increasing the size of the nuclear corrections. 
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